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ABSTRACT 

We report the observations of X-ray oscillations during the flare in a cool active star £ Boo for the 
first time. £ Boo was observed by EPIC/MOS of XMM-Newton satellite. The X-ray light curve is 
investigated with wavelet and periodogram analyses. Both analyses clearly show oscillations of the 
period of ~ 1019 s. We interpret these oscillations as a fundamental fast-kink mode of magnctoacoustic 
waves. 

Subject headings: stars: activity, - stars: coronae, - stars: flare, - MHD, - waves 



1. INTRODUCTION 

The magnetohydrodynamic (MHD) waves and oscilla- 
tions in the solar and probably in the stellar atmospheres 
are assumed to be generated by coupling of complex mag- 
netic field and plasma. These MHD waves and oscilla- 
tions are one of the important candidates for heating 
the solar/stellar coronae and accelerating the supersonic 
winds. In the Sun, magnetically structured flaring loops, 
anchored into the photosphere, exhibit various kinds of 
MHD oscillations (e.g. fast sausage, kink and slow acous- 
tic oscillations). The idea of exploiting such observed 
oscillations as a diagnostic tool for deducing physical pa- 
rameters of the coronal plasma was first suggested by 
Roberts et al. (1984) . Wave and oscillatory activity of 
the solar corona is now well observed with modern imag- 
ing and spectral instruments (e.g. instruments on-board 
/it SOHO, TRACE, HINODE and NoRH ) in the visible, 
EUV, X-ray and radio bands, and interpreted in terms 
of MHD wave theory (Nakariakov & Verwichte 2005) . 

Stellar MHD seismology is also a new developing tool 
to deduce the physical properties of the atmosphere of 
magnetically active stars, and it is based on the anal- 
ogy of the solar coronal seismology (Nakariakov 2007). 
The first stellar X-ray flare oscillations in AT Mic have 
been reported by Mitra-Kraev et al. (2005) and it was 
interpreted as a signature of slow magneto-acoustic oscil- 
lation. Mathioudakis et al. (2006) have reported the fast 
sausage oscillations of ~ 10 s in the flaring stars EQ Peg 
B. Mullan et al. (1992) have also found the oscillations 
on the scale of few minutes and interpreted as a transient 
oscillations in the stellar loops. In this Letter, we find 
the observational signature of the fast kink oscillations 
in the magnetically active star £ Boo (V = 4.55 mag). 
In the case of solar corona, such type of modes are im- 
portant to probe the local plasma conditions (e.g., van 
Doorsselaere et al. 2008). Therefore, such modes may 
also be useful to diagnose the stellar coronae. £ Boo is a 
nearby (6.0 pc) visual binary, comprising a primary G8 
dwarf and a secondary K4 dwarf with an orbital period of 
151 yr. In terms of the outer atmospheric emission, UV 
and X-ray observations show that the primary dominates 
entirely over the secondary (e.g. Laming & Drake 1999). 
The observations, data reduction, results and discussion 
along with the conclusions are given in the forthcoming 
sections. 



2. OBSERVATIONS AND DATA REDUCTION 

The star £ Boo was observed by RGS and EPIC-MOS 
on boards XMM-Newton satellite (Jansen et al. 2001) 
on 2001 January 19 at 19:25:06 UT and lasted for 59 ks. 
The detailed method of EPIC-MOS data reduction, and 
the timing, spectroscopic and flare analysis are given in 
Pandey & Singh (2008). Here, we used the 0.3-10.0 keV 
X-ray data observed from EPIC-MOS1. The X-ray light 
curve of £ Boo is shown in Figure [TJ where flare regions 
are represented by Fl and F2, and the quiescent state is 
represented by Q. The light curve was binned with 40s. 
In the X-ray light curve of £ Boo, a period of height- 
ened coronal emission 'U' was identified after the flare 
F2, where the average flux was 1.8 times more than the 
quiescent state. We have chosen the U region for our 
wavelet and periodogram analyses because of the stable 
evolution of the integrated 0.3-10 keV light curve after 
the flare F2. To study the emission in the U regions 
only, we fitted the spectral data using 1-tcmpcrature 
plasma model apec (Smith et al. 2001), with the qui- 
escent emission taken into account by including its best- 
fit 2-temperature model as a frozen background contri- 
bution. The best fit temperature and emission measure 
during the U part of the light curve are 8.7±0.3 MK and 
3.8 ± 0.5 x 10 51 cm~ 3 , respectively with reduced x 2 per 
degree of freedom of 0.99. During the flare F2, the av- 
erage values of temperature and emission measures were 
found to be 12.0 ± 0.6 MK and 3.5 ± 0.1 x 10 51 cm" 3 , 
respectively. 

Line fluxes and positions were measured using the 
XSPEC package by fitting simultaneously the RGS1 and 
RGS2 spectra with a sum of narrow Gaussian emission 
lines convolved with the response matrices of the RGS in- 
struments. The continuum emission was described using 
Bremstrahlung models at the temperatures of the plasma 
components inferred from the analysis of EPIC-MOS1 
data. The emission measure derived from the analysis of 
the EPIC data was used to freeze the continuum normal- 
ization. For the present purpose we give the line fluxes 
of He-like triplets from O vn in Table 1. 

3. ANALYSIS AND RESULTS 

3.1. Density Measurement 

The He-like transitions, consisting of the resonance line 
(r) Is 2 1 Sq - ls2p 1 pi, the intercombination line (i) 
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TABLE 1 

Observed O vii line fluxes during the Q and F2+U states 



Quiescent (Q) 


Flaring (F2- 


t-u) 


A(A) Flux Q 


A(A) 


Flux" 


22.08(f) 2.1 ±0.2 


22.08(f) 


2.2 ±0.3 


21.78(i) 0.8 ±0.2 


21.77(i) 


0.9 ±0.2 


21.58(r) 3.0 ±0.3 


21.57(r) 


3.4 ±0.3 



a Measured flux in 10 4 photons/cm 2 /s 
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FlG. 1. — The X-ray 40-s binned light curve of the star £ Boo in 
the energy band 0.3-10. keV. The overplotted continuous line is 
200 s binned light curve 
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Fig. 2. — The O VII line of £ Boo observed by RGS1. Left panel: 
Quiescent state 'Q' and right panel: Flaring state 'F2±U'. 



Is 2 % - ls.2p 3 P1, and the forbidden line (f) Is 2 % 
- Is. 2s 3 Si are density- and temperature-dependent (e.g. 
Gabriel & Jordan 1969). The intensity ratio G = (i+f)/r 
varies with temperature and the ratio R =f/i varies with 
electron density due to collisional coupling between the 
metastablc 2 3 S upper level of forbidden line and the 2 3 P 
upper level of the intercombination line. In the RGS 
wavelength, the O vii lines are clean, resolved and po- 
tentially suited to diagnose electron density and temper- 
ature (see Fig. 2). These lines at 21.60A(r), 21.80A(i) 
and 22.10A(f) have been used to obtain the temperature 
and density values from the G- and R-ratio using CHI- 
ANTI database (version 5.2.1; Landi et al. 2006). The 
G-ratio for £ Boo implies a temperature of 1.8 ± 0.6 and 
2.0 ± 0.6 MK during the quiescent and flare states, re- 
spectively. For these temperatures, we derive electron 
densities of n e = 1.3^g 2 x 10 10 cm~ 3 for the quiescent 
state (Q) and n a = 1.6^ $ x 10 10 cm" 3 for flare state 
(F2+U). The electron density was slightly higher during 
the flare state. However, this increase is well within the 
la level. 



3.2. Loop Length 

The length of the stellar loops associated with the ob- 
served flaring activity can be estimated from the various 
loop models. Below we derive loop length from four dif- 
ferent models. 

3.2.1. Hydrodynamic Flare model 

The details of this method are given in Pandey & Singh 
(2008). The derived loop length for the flare F2 is found 
to be of 7.9 ± 6 x 10 10 cm. However, the value of £ (=2; 
slope of density temperature diagram) is outside the do- 
main of the validity of the method. Therefore, the esti- 
mated loop length from this method may not be actual. 
Alternatively, Reale (2007) derived the loop length (L) 
from the rise phase and peak phase of the flare. The loop 
length of the flare F2 from this approach is estimated to 
be 3.6 ±0.8 x 10 10 cm. 

3.2.2. Pressure Balance method 

Shibata & Yokoyama (2002) assume that the gas pres- 
sure is balanced by the magnetic pressure and give 
the equation for the magnetic field (B) and the loop 
length(L): 

/ em \- 1/5 ( «<= v /l0 f t y 7/1 ° 

S_5 ° Vl0 48 cm--V V 10 9 cm- 3 ) \10 7 K ) (1) 

where n e is the density outside the flaring loop. By 
applying the average values of EM and T during the flare, 
we obtained a magnetic field strength of 29.4 ±6.0 G and 
a loop length of 3.6 ± 0.9 x 10 10 cm for the flare F2. 

3.2.3. Rise and decay time method 

The approach is described in Hawley et al. (1995). In 
this approach, the loop length has been derived by solv- 
ing the time-dependent energy equation in rise and decay 
phase separately in different limits: strong evaporation 
driven by heating in the rise phase and strong conden- 
sation driven by radiation in the decay phase. The time 
rate change of the loop pressure (p) for the spatial aver- 
age is given by 

\i> = Q-R (3) 

where Q is the volumetric flaring heating rate, and R is 
the optically thin cooling rate. The rise phase is domi- 
nated by strong evaporative heating, i.e., Q » R, while 
the decay phase is dominated by cooling and strong con- 
densation, i.e., R >> Q. At the loop top, there is an 
equilibrium, i.e., Q = R and the loop length can be given 



where is the decay time, r r is the rise time, Ta is 
the temperature at flare apex, and Xd = c q /c max . c max 
(=3.16 counts/s) and c q (=1.24 counts/s) are count rate 
at the peak of the flare and the count rate during the 
quiescent state, respectively. Using this approach, the 
loop length for the flare F2 is estimated to be 3.9 ± 0.5 x 
10 10 cm. The values of r d (=9885 s), r r (= 3351 s), and 
T A (=13.3 MK) are taken from Pandey & Singh (2008). 
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3.2.4. Haisch's Approach 

Given an estimate of two measured quantities, the 
emission measure (EM) and the decay time scale of the 
flare (ra ) Haisch's approach (Haisch f983) leads to the 
following expression for the loop length (L) 



fl/43/4 



(5) 



Using this approach, we obtained a loop length of 3.8 ± 
0.2 x 10 10 cm. The minimum magnetic field strength is 
estimated from the condition B 2 /8ir < 2n kBT, where 
fcfl is Boltzmann's constant. The minimum magnetic 
field strength thus obtained was 36 Gauss. This lies in 
the upper limit of the magnetic field determined by the 
pressure balance method. 

The loop length derived from the above all four meth- 
ods is consistent with each other. 

3.3. Wavelet and Periodogram Analyses 

We have used the wavelet analysis IDL code "Ran- 
domlet" developed by E. O'Shea. The program exe- 
cutes a randomization test (Linncll Ncmec & Nemec 
1985; O'Shea et al. 2001) which is an additional feature 
along with the standard wavelet analysis code (Torrence 
& Compo 1998) to examine the existence of statistically 
significant real periodicities in the time series data. The 
advantage of using this test is that it is distribution free 
or non parametric, i.e., it is not constrained by any spe- 
cific noise model such as Poisson, Gaussian, etc. Using 
this technique, many important results have been pub- 
lished by analyzing approximately evenly sampled data 
(e.g., Banerjee et al. 2001; O'Shea et al. 2001, 2007; 
Srivastava et al. 2008a,b; O'Shea & Doyle 2009; Gupta 
et al. 2009). 

The wavelet power transforms of the U part of the X- 
ray light curve of cadence 40 s (see the top-left panel of 
Fig. 3a) are shown in the middle left panel of Fig. 3a, 
where the darkest regions show the most enhanced os- 
cillatory power in the intensity wavelet spectrum. The 
cross-hatched areas are the cone of influence (COI), the 
region of the power spectrum where edge effects, due to 
the finite lengths of the time series, are likely to domi- 
nate. The maximum allowed period from COI, where the 
edge effect is more effective, is ~ 7363s. In our wavelet 
analysis, we only consider the power peaks and corre- 
sponding real periods (i.e. probability > 95%) below this 
threshold. The period with maximum power detectable 
outside the COI is ~2883 s with a probability of ~94% 
and the repetition of only ~ 3 cycles in the time series. 
However, another peak is visible at a period of ~1019 
s with a probability of ^99% and the repetition of ~11 
cycles. Therefore, only the periodicity of ~1019 s satis- 
fies the recently reported more strict criteria of O'Shea 
& Doyle (2009), i.e. at least four cycle of repetition over 
the life time of the oscillations. The middle-right panel 
in Fig. 3a shows the global wavelet power spectra of the 
time series from which the statistically significant period 
(~1019s) is selected. In Fig. 3a, the bottom panel shows 
the probabilities of the presence of two specific frequen- 
cies (or periods) corresponding to the first and second 
highest peaks shown in the middle-left panel as func- 
tions of the time after the start of the observations . The 
solid line represents the probability corresponding to the 
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Fig. 3. — (a) The wavelet result for U part of the X-ray light curve 
of £ Boo. The top panel shows the variation of intensity, the wavelet 
power spectrum is given in the middle panel, and the probability in 
the bottom panel, (b) Power spectrum from periodogram analysis. 



major power peak of time series data, while the dotted 
line corresponds to the secondary power peak. Wavelet 
analysis was also performed in the decay phase of the 
flare F2, and no significant periodicity was found. 

We have also performed the periodogram analysis 
(Scargle 1982) in the U part of the X-ray light curve 
(see Fig. 3b). In the power spectrum, the highest peak 
corresponds to a period of ^1005 s (probability ~ 75%), 
which is consistent with that obtained from wavelet anal- 
ysis. Therefore, only the period of ~ 1019 s appears to 
be statistically significant in our analyses. 

3.4. Fast-kink Oscillations 

The observed periodicity may be the likely signature 
of MHD oscillations associated with the stellar flare 
loops. The overdense magnetic loops are pressure bal- 
anced structures and may contain fast-kink and sausage 
oscillation modes with phase speed greater than the lo- 
cal Alfvenic speed (ui/k > va), and slow acoustic modes 
with sub- Alfvenic speed (co < va) under coronal condi- 
tions. 

The expressions for the oscillation periods of slow (P s ), 
fast-kink (Pfk), and fast-sausage (Pf s ) modes are given 
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Derived parameters 



TABLE 2 



Model 



Loop Length 
(10 10 ) cm 



Slow Mode 



Theoretically estimated Period 



Fast-Kink mode" 
00 



Fast Sausage Mode 

00 



Observationally 
Estimated 
period (s) 



Hydrodynamic 3.6 ± 0.8 1586 ± 353 1004 ± 391 313 ± 121 

Rise and decay 3.9 ± 0.5 1718 ± 222 1087 ± 374 339 ± 116 

Pressure balance 3.6 ± 0.9 1586 ± 377 1004 ± 407 313 ± 127 1019 

Haisch's approach 3.8 ± 0.2 1674 ± 92 1059 ± 343 330 ± 107 



a Large error bars are due to the large density error, b Loop width was determined by assuming a/L — 0.1 



as follows (Edwin & Roberts, 1983; Roberts et al., 1984; 
Aschwanden et al., 1999) : 

2L 2L f / Co \ 2 1 1/2 Lio 
Jc T jc I \v A J J VT 6 

Pfk = — =4tt 1/2 - ~ 205 *° V 9 s (7) 

jet j V.B 2 + B 2 / Bio 

D 27ra 3/2 / Po + p e \ 1 / 2 a 8v ^ 

Pfs = =47r' a — ?s 6.4 - s (8) 

f c k \Bl + Bl) Bio 

where L, a, p, and B are the loop length, the loop 
width, the ion mass density and magnetic field, respec- 

1/2 

tively. The ct = c va/{c 2 + v 2 A ) . where c D and va 
are the sound and Alfvenic speeds respectively. The Ck 
is the phase speed. Under coronal conditions, B « B e , 
p e << p , riH ~ n e , and c << va- The subscripts ' d 
and 'd stands for 'inside', and 'outside' of the loop. The 
j is the number of nodes in fast-kink and slow oscilla- 
tions, and is equal to 1 for the fundamental mode. The 
approximate relation in Equations 6-8 are given in coro- 
nal conditions for the fundamental oscillation periods of 
different modes in magnetic loops. Therefore, these ap- 
proximate equations are also valid for the stellar loops 
in the coronae of magnetically active Sun-like stars. The 
parameters are expressed in units of Lio = L/10 10 cm, 
a s = a/lifcm, n 9 = n o /10 9 cm- 3 , T = T e /10 6 K, and 
Bio = B/10G. The loop width was determined by as- 
suming a = 0.1L (Shimizu 1995) for a single loop. 

Using the observationally estimated loop length and 
width, density, magnetic field, and Equations (6)-(8), 
the estimated periods for different oscillation modes are 
given in Table 2. The observationally derived oscillation 
period ~ 1019 s approximately matches with the theo- 
retically derived oscillation period for a fundamental fast 
kink mode. 



such oscillations in the stellar loops. Recently, Cooper et 
al. (2003) have found that the transverse kink modes of 
fast magneto-acoustic waves may modulate the emission 
observable with imaging telescopes despite their incom- 
pressible nature. This is possible only when the axis of 
an oscillating loop is not orientated perpendicular to the 
observer's line of sight. 

Resonant absorption is the most efficient mechanism 
theorized for the damping of a kink mode in which energy 
of the mode is transferred to the localized Alfvenic oscil- 
lations of the inhomogeneous layers at the loop boundary 
(Lee & Roberts 1986). Intensity oscillations due to the 
fast-kink mode of magneto-acoustic waves do not show 
any decay of the amplitude during the initial 250 min- 
utes in the post-flare phase observations of £ Boo. This 
is probably due to an impulsive driver, which forces such 
kind of oscillations and dominates over the dissipative 
processes during the flare F2 of £ Boo. 

In conclusion, we report the first observational ev- 
idence of fundamental fast-kink oscillations in stellar 
loops during a post-flare phase of heightened emission 
in the cool active star £ Boo. These oscillations are a 
unique tool for the estimation of several key parameters 
of coronae, which are not directly measurable, e.g., the 
magnetic field. Some of the seismological techniques de- 
veloped for the solar corona, e.g., based upon the ratio of 
the periods of kink modes, do not require the spatial res- 
olution, and hence can be applied to estimate the density 
scale heights of the stellar coronae. Therefore, based on 
the solar analogy, such observations may shed new lights 
on the dynamics of the magnetically structured stellar 
coronae and its local plasma conditions. However, a fu- 
ture multiwavclcngth observational search is required to 
investigate various kinds of MHD modes to probe the 
local conditions of plasma in magnetically active stars. 



4. DISCUSSION AND CONCLUSIONS 

We found the first signature of fundamental kink oscil- 
lations with period of ~ 1019 s in the magnetically ac- 
tive star £ Boo observed by EPIC-MOS of XMM-Newton, 
which may be formed either by the superposition of op- 
positely propagating fast kink waves or due to flare gen- 
erated disturbances near the loop apex. The flaring ac- 
tivity may be the possible mechanism of the generation of 
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